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Abstract 

We solve exactly the scalar box integral using the Mellin-Barnes representation. 
Firstly we recognize the hyper geometric functions resumming the series coming from 
the scalar integrals, then we perform an analytic continuation before applying the 
Laurent expansion in e = (d — 4)/2 of the result. 



1 Introduction 



In 1974 't Hooft suggested that QCD in the planar limit was exactly solvable Unfor- 
tunately this program was too difficult to solve. There are however simpler quantum field 
theories in d — 4 that could enjoy this property. In particular N = 4 supersymmetric gauge 
theory is a cousin of QCD, ultraviolet finite, which is also a conformal field theory ( the one 
appearing in Maldacena's AdS/CFT correspondence j2]). Can we solve N = 4 super- Yang- 
Mills theory ? This is nowadays a hot topic, since its gravitational counterpart ( N = 8 
supergravity ) is also probably the only gravitational field theory which is ultraviolet finite. 

There are good reasons to believe that, in the 't Hooft planar limit, higher loop orders 
of N = 4 SYM are surprisingly simple. Indeed an iterative structure has been discovered at 
least in the four point planar amplitude, relating the two and three loop amplitudes to the 
one loop amplitude j3]-jl]-[5] and eventually allowing the perturbative series to be resummed 
into a simpler result. This result has been found thanks to a work of Smirnov [6], which 
has evaluated the Laurent expansion in e = (d — 4)/2 of the associated two- loop planar 
box integral. An important ingredient for his proof is the use of the Mellin-Barnes ( MB ) 
representation of the double planar box integral. However it is somehow disappointing that 
this technique works only perturbatively in the infrared regulator e. 

In this paper we analyze the MB representation in all cases where a non-perturbative ( 
in e ) amplitude can be computed, therefore concentrating our efforts to the four point one 
loop planar box integral, which is the basic information for the higher loops resummation. 
Also the arbitrary N(> 5)point one loop amplitudes can be reduced to a sum over a set of 
basic scalar box four point integrals. 

Back to the MB representation, the non-perturbative point of view we follow allows us 
to rewrite the harmonic series coming from the one loop amplitude in terms of generalized 
hypergeometric functions ( of several variables ) which in the massless and massive cases ( 
up to 2 external masses ) can be reduced to hypergeometric functions of a single variable. 

A key ingredient in our method is that it is always necessary to perform an analytic 
continuation of these hypergeometric functions before taking the Laurent expansion in e. 
We compare our findings with the direct method of integrating the scalar box and we find 
general agreement. 

Our non-perturbative method could in principle be generalizable to the study of the hard 
two mass scalar box integral [7] and eventually the double box scalar integral [6]- [8], where 
the hard part is to find the correct analytic continuation of an hypergeometric function 
of three and four variables ( respectively ) which cannot be reduced to an hypergeometric 
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function of a single variable. If we find the solution to this mathematical question we will 
report on it. 



2 Massless scalar box one loop integral 

Let us start with the massless scalar one loop box integral in d- dimensional space-time ( 
e = (d-2)/4: 



/ = T(2 — e) / dxi dx2 dx 3 dx± 8( x\ + X2 + £3 + £4 — 1 ) 
Jo 

[ xix 3 (-s) + x 2 x 4 (-t) Y' 2 (2.1) 

To proceed with the evaluation of the integral / we choose the Feynman parameters as 
follows: 



xi = (l-x)(l-y) 
x 2 = x(l-y) 

x 3 = yz 

x, = 3/(1 -z) (2.2) 
In this way the integral over y factorizes and we end up with: 



r 2 (e) 

r(2e) 



r(2 



dy (y(l-y)) 



1 f i 

dx dz [ (1 

n JO 
e-1 



x)z(-s) + (1 - z)x(-t) ] 



e-2 



r 2 (e) 

T(2e) 



(2.3) 



As we can see from eq. (I2.3p . the integration over x is straightforward and can be 
analitically performed. The resulting expression is: 



T 2 (e) f 1 

m T{1 - l) l dz 



z^i-sY" 1 



l-z 



{{l - z){-t) + zs) 



(2.4) 



By noting that the first integral over z stands for the Euler representation of the hyper- 
geometric function, we obtain the result: 
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dz z e - l {{-t) - z(-s - ty-^-sy- 1 - 1 ( S 



2 F X 1,6,6+1,1 + - (2.5) 

est V 1/ 



The other integral can be easily evaluated by replacing z — > (1 — z) 



dz z^H-s) - z{-s - t)Y\-ty- 1 = -t-^f 2 F 1 ( 1, e, e + 1, 1 + - 

6 S t \ S 



(2.6) 



By summing eqs. (12.51) and (12.61) we finally find that the massless scalar box integral is: 



r 2 (e)T(l- e ) 
r(2e) e s t 



-s) e 2 F! l,e,e + l,l + - + 



+ (-ty 2 F 1 1,6,6+1,1 + 



(2.7) 



The hypergeometric functions listed here admit the following Laurent expansion in the 
infrared regulator e: 



,l\ ( 1, e, e + 1, 1 + |) = 1 - e log (-0 - e 2 Lz 2 (l + I) + 0(e 3 ) 



producing the well known perturbative result 



(2.8) 



2 r(l-e)r 2 (l + e) 



ste 2 T(l + 2e) 
that can be ultimately simplified by noting that 



-sy +(-ty + e 2 Li 2 (--) + e 2 Li 



,7T 



— 6 



(2.9) 



-log 2 (-)-- 
2 B \tJ 6 



(2.10) 



There is a second non-perturbative expression, owing to the following identity 



- 2 F 1 (1, e,e + l,z) = - + — !— z 2^1 (1, 1 + e, 2 + e, z) 
e e 1 + e 



(2.11) 



which allows to perform more easily the perturbative expansion in e 
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r 2 ( e ) frn 



+ 
+ 



r(2c) 
r(i - *) 
(i + 6) 



s t 



s t 



+ 



s t 



1 + 2 F 1 (l,l + e,2 + e,l + ^ 



s t 



1 + -M 2 fi(l,l + e,2 + e,l + ^ 



(2.12) 



Also this second hypergeometric function can be developed in power series of e: 



1 + e 



Z 2^1(1,1 + 6,2 + 6,2) 



- loef 1 - z) 



eLi 2 {z) +0(e 2 



(2.13) 



In this case we note that in the second part of eq. (I2.12p the term of order zero in e 
cancels out while the term of order e is given by 



Lz 2 (l + - Lz 2 fl + ~) - log(-s)log(-^) - log(-t)logf-^ 
Li 2 (- S -) + Li 2 (--) - (2.14) 



tJ V s , 

In this way the formula ( 12. 9 \ is recovered. 



3 Mellin - Barnes representation of the massless scalar 
box integral 

We are going to compare the non-perturbative calculation of the massless scalar box (12. 7p 
with the MB representation, a powerful method that has allowed to perform analytically 
more complicated integrals, like the scalar double box integral. 

The MB representation allows us to rewrite the massless scalar box integral 



r 2 (e) 



r(2e) 
in the following way 



r(2-e) f dx f dz [ (1 - x)z (s) + (1 - z)x (-t) ] e ~ 2 (3.1) 
Jo Jo 
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T(2 - e) [ (1 - x)z (s) + (1 - z)x (-f) Y~ 2 = 

= — dw T(2 - 6 + w) F(-w) , , ( (1 f (3.2) 

2m 7_ JOO 1 ; 1 J ((i_ x ) z (_ s ))2-^ ^ > 

The integrals in x and z are now factorized and produce: 



/ dx x w (1 - x) 
Jo 

f 

Jo 



e-2-w 



F(w + l)F(e - 1 - w) 
.,_ 2 _„. Kir I 1)1 U 1 w) 



dz{l-zYz^- w = r( v e) L (3.3) 

from which we obtain the MB representation of the massless scalar box as 



1 f +io ° dw (-t) w 
I = ^7 , { J-^„ r > + 1) T(2 - e + w) F(-w) r 2 (e - 1 - w) (3.4) 



r(2e) J_ ioo 2m (-s)^+ w 

We can close the integration contour to the left or to the right and the result is unaffected 
by this choice. We always elaborate these integrals exactly without developing in Laurent 
expansion of e; our point of view is therefore different from the usual applications of the MB 
representation we have seen in literature. 

For example, let us choose to analyze the poles to the left: 

1) it is more simple to analyze the simple poles w = e — 2 — n; the corresponding series 
can be resummed into the following hypergeometric function 

/ - r2(e) r2(1 ~ e) ( tY~* F (l 1 2 e S ) (3 5) 

At fist sight this result is not much interesting, since this hypergeometric function doesn't 
admit a meaningful development in e like those of the direct method. It turns out that it is 
necessary to perform an analytic continuation into the right variables, before that developing 
in e makes sense. 

We are forced to choose the following analytic continuation 



> F > = H) r ( i + ort-i) ' F ' + + i 



sj r(e-l) V t 
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It then appears the first part of the standard solution obtained with the direct method, 
( see eq. flU}) 



I (2e) est \ s J 

and a spurious term 

r 3 (e)r 2 (l-e) f ( s\-* 

' s e (1 + -) (3.8) 



T(2e) st V t 

2) let us consider now the contribution of the double poles. It turns out that in order 
to get a well defined resummed expression it is necessary to regularize the double poles 
introducing a parameter 5: 



l r +i °° dw (—t) w 

1 = m&L ^(^ ;r( '° +1+ ' )r( '° +1)r(2 - £+ro)r( - m) 

r(e-l-w) T{e-l-w-8) (3.9) 

Case 2a) first consider the simple poles w = — 1— n — 5. Their contribution can be easily 
resummed into the following formula 



A development in 5 is necessary 



T(e + 5)T(l-e-5) = T(e)r(l - e) (1 + 5(^(e) - i>(l - e))) + 0(6 2 ) (3.11) 
It follows that 

s^o 5 r(2e) V st \ t) 

- W) T(l -' ] KJ-^-^+^f) *=r ( 1 + 7)"' + °w < 3 - 12 > 

12a contains a singular term in S and a finite term, which is the main contribution. 
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Case 2b); we analyze now the simple poles w = — 1 — n. Their contribution can be 
summarized as 



s "> 



T 2 (e) i r(e - 5)V(1 - e) (- 8 )« / 

hb = F(2l) & 5 r(i- W e) "7T" 2Fl I 1 ' e - *' 1 - 5 ' - 1 J ^ 3 - 13 ) 

Again we must choose a convenient analytic continuation in the variable (1 



(3.14) 



The firs term goes like 5 and in this case it is possible to take directly the limit (5—7-0: 

^-Si^^^(w + M + f) ,3,5) 

representing the residual part of the exact solution ( eq. (12.71) ). 
The second term can be resummed as 

lim^ra-e) t£l! ( 1 + £r e (3.16) 
5->o 5 r(2e) V ; s t V t/ V t/ v ; 

that can be decomposed in a term divergent in 5 plus a finite term 



Italia - e )fc^! fi + f 

<5^o 5T(2e) v s t V £, 

Collecting our partial findings, we obtain with MB representation the exact solution plus 
the following spurious terms 



r 2 (e) Wl (s) 

r(2e) 



r(l-e)i-^- + {log(-l) + r(c)r(l - e) ( cos(7re) - (-1) £ )} = 

(3.18) 



Surprisingly this lengthly expression is zero due to non trivial cancellation between all 
terms. 
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4 The scalar box integral with one external mass 



The next step is comparing the MB representation with the direct method in the case of 
one external mass. To be precise we start with the following integral 



I m (s, t, m) = T(2 — e) / dx\ dxi dx 3 dx± S( X\ + X2 + £3 + X4 — 1 ) 

Jo 

[ xix 3 (— s) + x 2 x 4 (— t) + Xix 4 (— m 2 ) ] e 2 



(4.1) 



We choose the Feynman parameters as in the formula( l2.2p and we find the following 
expression 



1 rl 



I™ = Ii| r(2-c) J dx j dz[{l-x)z {-a) + (1 - z) x (-t) + (1 - x)(l - z)(-m 2 ) p 2 



(4.2) 



The integration in x is again obvious and it gives rise to 



T 2 (e) f 1 



U-j) + (1 )-'-((!- z)(- t)) 

(1 — 2;)(m 2 —t) + zs 



\e-l 



(4.3) 



It is convenient introducing the following notations 



z 



m 2 — t 
m 2 — t — s 



m 



z\ 



m — s 



from which the integral I m we want to compute is of the form 



(4.4) 



r 2 (e) 

r(2e) 



T(l-e) 



,m — s) 



-1 /-i 



dz 



z - Zi 



ie-1 



(s + t - m 2 ) y 2; - z 



:-t) 



e-l /•! 



,e-l 



(s + 1 — m 2 ) J z - z 



(4.5) 



The first part can be rearranged as follows 



^ dz {Z ~ Z%r ~ X = {- Zl ) 



e-l 



z- Z 



Z\ 



Z\ — Zq J J 



dw w 6 ' 1 I 1 



Z\ 



-XV 



Z\ — Zq 
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e-1 



€-1 



1-^1 

<2o — ^1 
«1 



dw w*- 1 I 1 - 



-1 



(l-^i) 6 - 1 fl-zt 



Zq — Z\ 



Fx l,e,e + l, 



Zq — Z\ 
Zl 



Fx 1, e, e + 1 



Z\ — Zq 
1- Zl 



Zq — Z\ 



(4.6) 



By substituting eq. (14.61) in I m and taking into account the definitions of Zq and Z\ in 
terms of m 2 , s, t we finally obtain 



rrrj 
1 l 



r 2 ( e ) r(i-e) f (-5 



T(2e) 6 



Fi ( l,e,e + l, 



s + t — m z 



—m 

s t 



Fi l,e,e+l 



m 2 (s + 1 — m 2 ) 
si 



(4.7) 



The second part of eq. (14. 5p is easily computed by replacing z — >■ (1 — z) 



r 2 (6) r(i-e) 

r(2c) e 



s t 



*) e / s + t-m* 

2 Fi l,e,e + 1, 



The complete result is simply the sum of eqs. fl4.7p and (14. 81) 



T m = jm + I m 



(4.9) 



In the limit m? — > the massless scalar box integral I(s,t) is found. 

The result (]4.9[) can be put in another form which makes more explicit the Laurent 
expansion in e of the non-perturbative solution 



+ 



r 2 (e) fr(i-e) 

r(2e) 

r(i- 



-777. 



2\e 



S t 



S t 



S t 



(1 



— sY ( s +t — m z 



s t 



t 



'-tY ( s + t-m 2 , , , 

+ M~ f I 2^(1,1 + ^2 + 6, 



s t 



'—m 2 Y {m 2 (s + t — m 2 ) 



2 F X ( 1, 1 + e, 2 + e, 
s + t — m 2 



s + t — m z 



s t 



s t 



Fx 1,1 + 6,2 + 6 



m 2 (s + t — m 2 ) 
si 



(4.10) 



In the expansion in e it is enough to recall that 
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1 z 2 F 1 (l,l + e,2 + e,z) = -log(l-z) - e Li 2 {z) + 0{e 2 ) (4.11) 



1 + e 



The logarithmic terms cancels out between them, and at the order e we find using the 
identity 



TV 2 



Li 2 (l-x) = Li 2 (x) - — + log(x) log(l-z) (4.12) 

o 



the following contribution 



,'m 2 — t\ r { m 2 — s\ ( (m 2 — s) (m 2 — t)\ it 2 
Li 2 + Li 2 ( — - — - Li 2 ( i 1 - — (4.13) 



st J 6 

5 Comparison with the Mellin - Barnes integral repre- 
sentation 

The integral I m (eq.f l4.ip ) contains the sum of three terms and can be developed using the 
MB representation by decomposing twice that sum 

r(2 — e) ( x-lXs{— s) + x 2 x 4 (-t) + Xix 4 (-m 2 ) ) e ~ 2 = 

(x^i-m 2 ))* ( Xl x 3 (-s)y (x 2 x 4 (-t)y- 2 - a -P (5.1) 
Introducing the usual Feynman parameters ( see eq. (12.21) ) we find 

i r +i °° do r +i °° dB 

/m = / aw / =£:r(-a)r(-/3)r(2-e + a + /3)r(e-l-a-/3) 

T(l + p) T(e - 1 - /3) T(l + a + /3) (-m 2 ) a (-s)^ (-t) e - 2 - a -P (5.2) 

It is more convenient integrating firstly in a and then in /3 

{-ty- 2 dp , 8 \p 



r(2 



/•+*°° /ex/3 

L ^ (?) r (-« r ( i +« r (- i -« 
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r +iQO da fm 2 \ a 

/ — — r(-a) r(e-l-a-/3) T(2-e + a + f3) T{1 + a + (3) 

J-ioo 2vn \ t J 

(5.3) 

1) let us discuss the poles to the left of a 

a = -l-k-p (5.4) 
The corresponding series can be resummed as 

Without loss of generality, we close the contour integral in f3 to the right. There are two 
contributions 

Case la: (5 = n 

-P^-iF, (l,l,2-e,^L-) (5.6) 



1 — e \ m 2 — t 



and 

Case lb: /3 = e — 1 + n 



m 2 — t 



I -+ T 2 (e)T(l-e) ( — ] ( _ J (5.7) 



m 2 — s — t 



By using the formula of analytic continuation 



r 2 (6)r(i - e) _ 



r(c-i) 

the whole contribution simplifies to 



z 6 - 1 (1 - z)~ e (5.8) 
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therefore the case 1, which is the simplest one, results in 



h " T^f(27) 1 V ' ' ' s ) (5 ' 10) 

2) let us analyze the more difficult case, i.e. the series of poles 

a = -2-k-/3 + e (5.11) 
r(2e) ^ 1 ^ 1 ! WJ 

r°° ^r- (-^Y r (-/3) r(i + P) r(e - 1 - p) v{k + 2 + p - c ) (5.12) 

Again we choose for simplicity to analyze the poles to the right 

Case 2a: (3 = n 



ha ~ r(2e) MC M f )L r(H2 - £ )r(n + 2- f ) W \w?) { ] 



This double series is tabulated jS] and corresponds to a generalized hypergeometric func- 
tion of two variables ( in particular of type 2) 



= _ r(or(i - or(e-i) Fj L _ 2 _ 2 _ t . , 

1 (2eJ \ rrr m z 



Let us note that this particular hypergeometric function of two variables can be reduced 
to a hypergeometric function of a single variable [9] 



F 2 (a,P,P',a,a;x,y) = (1 - x)^ (1 - y)^' 2 F X 1^3,(3', a, _ * V _ - ) (5.15) 
from which we obtain 
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rm (-™ 2 Y r(£)r(i-6)r(e-i) / St 

1 2a — ~ 7 — 5 Iw — 5 \ 777T\ 2-ri 1, 1, z - e, 



m 2 — t) (m 2 — s) T(2e) V ' ' ' ( m2 ~ s ) ( 777,2 ~~ t) J 

(5.16) 

It is necessary applying again the analytic continuation (|5.8p . to obtain 
(-m 2 ) e ^ 2 ' 



rm 
1 2a 



S t 



m^> 2Fl ( , i, e , e+ i,'"'(' + *-"■•)') 

r(2e) e V st / 



(-m 2 )- r 3 ( £ )r 2 (l-e) ^ m 2 (m 2 - g -t) V £ 
+ "Tf f(2e) I 7* ) (5 ' 17) 

This formula contains a part of the exact solution plus a spurious term, as in the massless 
case. 

Case 2b) let us analyze the poles (3 = n + e — 1. 

In this case the resummation is easier and results in 

IS = _ M_ iWi^ , Fl ( lf !, 2 _ <, _«_) ,5,8) 
s (m z — s) 1 (2e) (1 — e) \ m z — s J 

Applying again the analytic continuation formula ( eq. fl5.8p ). we obtain 



- s y r 2 ( g ) r(i-6) / s + t-m 2 
JT r(2e) e 2Fl ^' e ' e + 1 ' 1 

-s) £ r 3 (e)r 2 (:i - e) fm 2 -s-t 



s t r(2e) V * 

To summarize, considering the sum of the three terms 



(5.19) 



i m = IT + IZ + J£ (5.20) 

the spurious terms cancel out and we faithfully reproduce the exact result ( see eq. (14. 9jl ). 

To conclude the Mellin-Barnes method is completely equivalent to the direct method in 
the massive case, as in the massless case. We expect the same conclusion also in the case 
of easy two mass scalar box integral [7J, whose exact result can be represented in terms of 
hypergeometric functions depending only on a single variable. 
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6 Conclusions 



Using the Feynman parameter method, we have compared the Mellin-Barnes representation 
and the direct method of solving one loop scalar integrals at a non-perturbative level in the 
infrared regulator e. Usually the MB representation is used in literature perturbatively in 
e; in our approach we firstly recognize the hypergeometric functions resumming the series 
coming from the scalar integrals, then we perform an analytic continuation in the right 
variables before applying the Laurent expansion in e of the result. 

Our method has been tested in the massless and massive ( with one external mass ) cases, 
but in principle it could be generalizable to more complicated integrals, like the hard two 
mass scalar box [7J and the double scalar box [6|- [8]. Work is in progress in this direction. 
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